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Abstract 

We present an analysis of isospin-breaking effects in threshold pion-nucleon scattering due 
to the light quark mass difference and the dominant virtual photon effects. We discuss 
the deviation from various relations, which are exact in the isospin limit. The size of the 
isospin-violating effects in the relations involving the isovector ttN amplitudes is typically of 
the order of one percent. We also find a new remarkably large effect (~ 40%) in an isoscalar 
triangle relation connecting the charged and neutral pion scattering off protons. 
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1. Pion-nucleon (ttN) scattering is one of the prime reactions to test our understanding of 
the spontaneous and explicit chiral symmetry breaking QCD is supposed to undergo. During 
the last years, there has been considerable interest in using irN scattering data to extract 
information about the violation of isospin symmetry of the strong interactions [1, 2], some 
analyses indicating effects as large as 7% [3, 4]. In both these analyses, the source of this 
rather large effect remains mysterious. Microscopically, there are two competing sources of 
isospin violation, which are generally of the same size, namely the strong effect due to the 
light quark mass difference rrid — m u ~ m u and the electromagnetic one caused by virtual 
photons. For neutral pion scattering off nucleons, these effects can be dramatically enhanced 
due to the smallness of the isoscalar pion-nucleon amplitude [5, 6]. This spectacular effect in the 
difference of the 7r°p and 7r°n scattering lengths is, however, at present not amenable to a direct 
experimental verification. It is therefore mandatory to include also the channels with charged 
pions in any analysis of isospin violation. To do this in a consistent fashion, one has to develop 
an effective field theory (EFT) of pions, nucleons and virtual photons. The corresponding 
effective Lagrangian was developed in refs.[7, 6] extending the standard irN EFT (for a review, 
see ref. [8]). The pertinent power counting of the EFT is based on the observation that besides 
the pion mass and momenta, the electric charge e should be counted as an additional small 
parameter, given the fact that e 2 /47r ~ M 2 /(4ttF 7T ) 2 ~ 1/100 (with M n and F w the pion mass 
and decay constant, respectively). From here on, we collectively denote these small parameters 
by q. Similar information can also be obtained from precise data on pion photoproduction, as 
detailed in ref. [10] (for an overview, see the talks [11, 12]). The aim of this paper is to give a 
first systematic study of the expected size of isospin violation in the ttN amplitude at threshold 
based on a set of relations, which are fulfilled in the limit of exact isospin. We stress again that 
in the framework we are using, a consistent separation of the electromagnetic and the strong 
effects is possible and to our knowledge this has not been achieved before. Only when a mapping 
of the method developed here on the commonly used procedures of separating electromagnetic 
and hadronic mass effects (such as the NORDITA method [13]) has been performed, a sensible 
comparison with the numbers quoted in the literature will be possible. 

2. Consider now the process 7r a (g a ) + N(p±) — > 7r 6 (cfo) + N{jp2), where 7r a denotes a pion 
of (cartesian) isospin a and N the nucleon. In the centre-of-mass system (cms), the four- 
momentum of the incoming nucleon is m • v + p\ = (E\ = m + v ■ p±, —q* a ), the one of the 
outgoing nucleon is m-v + p2 = (E2 = m + v -p2, — (ft), where m denotes the nucleon mass in the 
chiral limit. Similarly, the incoming pion has q a = {uj a ,q a ) and the outgoing pion q^ = (101,,%). 
The analysis of isospin violation in irN scattering proceeds essentially in three steps. First, 
one ignores all isospin breaking effects, i.e. one sets e = and m u = m^. Only if within this 
approximation one is able to describe the low ttN partial waves in the threshold region as given 
by various partial wave analyses, one can be confident to have a sufficiently accurate starting 
point.* 5 That this is indeed the case was demonstrated in refs.[9, 14, 15]. Ref. [15] also contains 
a detailed discussion of the kinematics pertinent to the case considered here. In the second step, 
one should include the leading isospin breaking terms encoded in the pion and nucleon mass 
differences. The corresponding terms in the effective Lagrangian read [7, 6] (also shown are the 
terms responsible for the explicit chiral symmetry breaking) 

C*n = d^dV + B (m u + m d )7r 2 + 2Ce 2 TT + 7r- + ... , 

#5 Note that in the available partial-wave analyses, electromagnetic and some hadronic mass effects are generally 
removed by some methods like e.g. the one from NORDITA [13]. 
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C nN = N (ido + ci M 2 (l + vr 2 /F 2 ) + B c 5 (m u - m d ) r 3 (l + vr 2 /F 2 ) + / 2 e 2 r 3 + . . .) AT,(1) 

where Bq is related to the scalar quark condensate, Bo = \ (0\qq\0)\/F%, and F is the pion decay 
constant in the chiral limit, F w = F[l + 0(q 2 )], The pion mass difference M 2 + — M 2 is entirely 
determined by the low-energy constant (LEC) C [16], * 6 while to third order in small momenta 
the strong (electromagnetic) proton-neutron mass splitting is given by the LECs C5 (fz) [6]. Note 
that the operator ~ C5 does not only contribute to the strong np mass splitting but also has a 
contribution ~ NNirir to the two-pion vertex which will be of relevance later. In terms of the 
operators defined in eq.(l), retaining only the terms leading to the strong and electromagnetic 
(em) hadron mass splitting is achieved by setting 

Ce 2 / , f 2 e 2 + , m u - m d / , but e 2 = . (2) 

This is the approximation which we will consider here. In fact, in neutral pion photoproduction 
off nucleons, to third order in small momenta, this approximation leads to the only isospin 
breaking effect, which reveals itself in the large cusp effect at the secondary threshold (i.e. at 
the 7r + n threshold in the case of 7^ — > 7r°p). In the third step, which goes beyond the scope of 
this paper, one has to account for all virtual photon effects, in particular soft photon emission 
from charged particle legs and the Coulomb poles due to the ladder exchange of (hard) virtual 
photons between charged external particles. In that case, the notion of partial waves becomes 
doubtful and one better compares directly to the available cross section and polarization data. 
We believe, however, that the essential effects of isospin violation are captured in the calculation 
presented here. 

3. In the presence of isospin violation, i.e. isovector symmetry breaking terms such as {m u — 
md) (uu — dd) , one has to generalize the standard form of the ttN scattering amplitude to 

T ab (u,t) = rt+( W> t) + ^t*^ , (3) 

in terms of two isoscalar (T ah ' ) and two isovector amplitudes (T ab ' ). These are functions of 
two variables, here we choose the pion energy oj and the invariant momentum transfer squared 
t. More precisely, oj can be chosen to be either the energy of the in-coming or out-going pion, 
since these are no longer equal 



Auj = UJ b - 0J a = 



(M 2 - Ml) - (m 2 - m 2 ) _ (M 2 - M 2 ) - (m 2 - m 2 ) r 



mi 



2m\ 

(4) 

with M a b (mi^) the mass of the in-coming, out-going pion (nucleon) and ^/s the total cms 
energy. It is important to note that while the pion energies iv a ^ are of order q, their difference 
only starts out at second order in the chiral expansion. This has important consequences as will 
be discussed later. The T-amplitudes split, of course, into a spin non-flip and a spin-flip term, 
denoted by g and h, respectively (for more precise definitions, see e.g. ref.[15]). At threshold, 
only the spin non-flip amplitudes can contribute and eq.(3) simplifies to 

T? N , thl = { 5 ab g+ + 6 ab r 3 git + i e bac r c g~ b + i e bac r c r 3 <fc } , (5) 



#6 Note that we work in the cr-model gauge so that the term C{QUQU') only contributes to two-point functions. 
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with Ni = y/{Ei + m i )/2m i (i = 1,2) the standard spinor normalization (its relevance is dis- 
cussed in detail in ref.[17]). In what follows, we do not consider these normalization factors 
since they are related to the external kinematics. Consequently, the isospin violating effects are 
entirely confined to happen within the given Feynman graphs we consider. The g a £ are, of 
course, proportional to the corresponding S-wave scattering lengths. Note that in the presence 
of isospin violation, these amplitudes can become complex even at threshold (since the mass of 
the incoming two-particle system is no longer equal to the mass of the out-going one). It is also 
important to realize that the g^ b are exclusively sensitive to the neutral to charged pion mass 
difference (i.e. the LEC C) whereas the are given by the operators ~ f 2 , ~ C5, i.e the strong 
and em proton-neutron mass difference. 

4. Isospin violation is best characterized in terms of quantities which are exactly zero in the 
isospin limit of equal quark masses and vanishing em coupling. With the three pion (tt^, it ) and 
two nucleon (p, n) fields, we have ten reaction channels, which in the case of isospin symmetry are 
entirely described in terms of two amplitudes. One thus can write down eight isospin relations 
(see also ref.[18] for a general analysis) 



#1 = 2 



T^+p — >7r +p ~\~ T^—p — > n —p ~\~ 2 T^Op^^Op 

2 9li + 5 2 + 2 - 2 533 + g ft + oli - 2 olt (6) 
ali + al 2 + 2533 + aft + all + ^all ' 

j-, n -^7r+p— »7r+p T n -p^ n -p — \^2T 7T -p^ 7T o n 

1X2 = £ 7= 

-^7r+p— »7r+p -^7r — p— >n~p V 2T 7r -p^ 7r o J! 

9i2 + 921 ~ 9i3 ~ 323 + 9u + all - 9i3 - 323 
9l2 + 921 + 9l3 + 523 + 9l2 + dll + 5?3~ + 323 



(7) 



R 3 = 2 



rp rp 



_ 2 531 + 532 5_13 523 + 531 + 532 5_13 523 ^ 
531 + 532 + 313 + 523 + dll + dl2 + 013 + 52s" ' 

R4 = 2 ^r^P ~ T ^-n = 2 dtl + + 3ll ~ 3I2 ~ 5jT (9) 

J 7r+p^7r+p ~r J 7r _ n^7r-ri 511+522 5l2 521 

3+ _i_ 3+ I 3— _i_ 3— 

R 5 = 2 ZlZIZHZl 1 TT+n-^Tr+n _ 2 5n + 522 t 5l2 t 521 

^7r _ p— »7r _ p "I" -^7r+n-+7r+n 5li + 522 + 5l2 + 521 



T -1 T -1 3+ 
1 7T°p— >7T°p _ 7T^7l — >7T^TL 2 ^33 

^7T°p— >7T°p ~l~ -^~7r rt^7r°ri 533 



^ _ 2 Tt u p^n u p j 7r u n^7r u n _ 2 ^33 



^7 = 2 * p ^ n " = fl 3 , (12) 

^7r p— +7r (, n ~r ^ 7r+n— >7r°p 

D r, - 1 7T°p— ^T" 1 "" J 7r°?l— >7T-p D 

^8 - ^ ™ —7^ — = --K3 • {i-O) 

1 TT°p-^1T + n ' 1 TT n^W p 

The first two, the so-called triangle relations, are based on the observation that in the isospin 
conserving case, the elastic scattering channels involving charged pions are trivially linked to the 
corresponding neutral pion elastic scattering or the corresponding charge-exchange amplitude. 
To be precise, these ratios are to be formed with the real parts of the corresponding amplitudes 



3 



evaluated at the pertinent threshold, symbolically T %a ^^ % b^ should read Re T^ N ^b N - The 
imaginary parts of some of the amplitudes will be discussed later. Of particular interest is 
the second ratio, which is often referred to as the triangle relation. Only in this case all three 
channels have been measured (for pion kinetic energies as low as 30 MeV in the cm system) 
and the 7% strong isospin violation reported in refs.[3, 4] refers to this ratio. We stress again 
that it is difficult to compare this number with the one obtained in our calculation since a very 
different method of separating the em effects is used. The ratio Rq parametrizes the large isospin 
violation effect for 7r° scattering off nucleons first found by Weinberg [5] and sharpened in ref.[6], 
Rq ~ 25%. Note that in Ri the isovector terms drop out completely and one thus expects also 
a large isospin violation in this ratio (since the isoscalar parts are strongly suppressed and are 
of the same size as the symmetry breaking terms). To our knowledge, this is the first time that 
this particular ratio has been called attention to. From an experimental point of view, it has 
the advantage of avoiding the almost unmeasurable nir° amplitude appearing in Rq. However, 
both i?i and Rq are sensitive to the precise values of the combination of LECs C2 + C3 — 2c\ 
since the strong contribution to the isoscalar scattering length is not even known in sign at 
present. The predictions for the other ratios are more stable since they involve the larger (and 
better determined) isovector quantities. Note that the relations Rj = R3 = — R$ follow from 
time reversal invariance. In what follows, we will calculate the six ratios Ri to leading one loop 
accuracy, i.e. to third order in small momenta. For that, we have to consider tree graphs, some 
with fixed coefficients and some with LECs, and the leading one loop graphs involving lowest 
order couplings only. 

5. The pertinent Born graphs calculated to first, second and third order are depicted in fig.2 
of ref.[15]. The ones contributing here are lb, 2a, 2b, 2d, 3a-3f and the additional tree graphs 
~ /2>~ C5 are shown in fig.l. Before giving the results for the sum of all Born graphs (tree 
graphs with or without LECs), some important remarks concerning the chiral power counting 
are in order. Although the so-called Weinberg-Tomozawa NNirir contact graph gives a first 
order contribution to g~ b , in the ratios Ri its effect is always proportional to Aw, which is of 
second order, cf. eq.(4). Consequently, isospin violation only starts at second order in the chiral 
expansion. Furthermore, in some tree graphs with intermediate nucleon lines the pion energy 
difference enters since v ■ p\ — v ■ P2 = ujb — uj a , which is thus also of order q 2 . We have only 
accounted for this difference whenever it was necessary and consequently neglected it when it 
would only lead to a fourth (or higher) order contribution. This applies in particular also to the 
loop graphs given below. Therefore, the final results depend on the choice of taking either the 
energy of the in-coming or the out-going pion as reference energy. This difference is, however, 
beyond the accuracy we are working and thus gives us the possibility to estimate some higher 
order effects. After mass and coupling constant renormalization, the sum of all Born graphs 
gives the following contributions (for notations, see ref.[15]) 

F$g+ = -4 Cl M 2 + 2(c 2 -|i)u, 2 + 2 C 3W 2 , (14) 
F 2 g s + = -2B (m u -m d )c 5 5 a3 + ±e 2 f 2 F 2 (5 a3 -l) , (15) 

uj{l-4^f 2 2 \ g^uj 3 u 1 x/,^|2 , , 1r n 
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= {-Bo{m u -m d )c b + \e 2 f 2 F^ - S b3 ) 

+ ^(Wa\ 2 + |%| 2 ){2 J B (m u - m d )c 5 + ^e 2 f 2 F^(5 a3 - 5 b3 ) . (17) 
At third order, we have to consider the set of one loop graphs shown in fig.2. These give 
FU + ab = -^(J^) + J c {-u>)){5 ac -5 cc )- 9 ^M%(J^ 

- & (m 2 (2M 2 - M 2 ) 7o -(0) + (M 2 - 3M 2 ) Jg(0))) - ^M 2 J«(0) , (18) 
FU~ ah = ^(4^(J C H- J C (- W ))+6A^) 

+t{K + + {Ml - M 2 )(A« - At) + t2 " 3t ^ 2 +M&2) ) 

„2 r (M 2 - m 2 ) 2 

+ H |4c[5t - 4(M a 2 + M 2 ) - 1 s - fej ]/ afc (t) + 20u,(At + A*) 



T (M 2 - M 2 )(A« - A^) - —(t - 3(M 2 + M 2 )) 



+ 



| (6 M Ai + At + At) - ^ (M 2 J d '(0) - At) (25 dc - <5 d 

+ ((f f\(M 2 Jm - At) + 8M 2 d 28 (A) - jrf c ^ , (19) 

where one has to generalize the standard loop functions Jo, A,,-, 70, Jo (see the appendices in 
refs.[8, 19]) to the unequal mass case as indicated by the superscripts 'a,b,c' referring to the 
pion isospin indices. The loop contributions to g~ b are divergent. These divergences are cancelled 
by the appropriate dimension three operators from C^ cm as constructed in ref.[6]. Note that 
in the approximation we are using, the finite parts of these terms are set to zero. We remark 
that the imaginary parts generated by the loops are very small for the threshold kinematics 
considered here. This is made more precise in the following. 

6. We are now in the position to analyze the ratios Ri as defined in eqs.(6-13). We use the 
standard masses as given in the PDG tables for m p . n , M w o,± [20]. For the mass parameter m we 
can use m = (m p + m n )/2 to the order we are working and we identify the leading term in the 
quark mass expansion of the pion mass with the neutral pion mass, Mq = M w o. We also use 
F w = 92.4 MeV #7 and g A = 1-26. The LECs ci j2 ,3,4 and <k are taken from fit l,2,3 #s of ref.[15] 
(we refer to that paper as "FMS"), C5 from ref.[9] and fa from ref.[6]. 

In table 1, we give the results for the ratios Ri that are not entirely given by isoscalar 
quantities. These numbers should be more stable than the ones for the isoscalar ratios -Ri,6- 

#7 In principle, we would have to differentiate between F n + and F n o. This difference is of second order and 
would therefore show up as a third order contribution due to the Weinberg-Tomozawa term. At present, the 
empirical determinations do not allow to differentiate between these two values and we thus work with one value 
given by the charged decay constant. 

#8 Fit 1,2 and 3 refers to the Karlsruhe [21], the Matsinos [22] and the VPI [23] partial wave analysis, in order. 
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This is indeed the case, there are no large variations between the three parameter sets given 
in FMS. For the pion energy u, we have used the arithmetic mean of the in-coming and out- 
going energies (this is also the most natural choice since it preserves the time reversal invariance 
between ir + n <-> ir°p and ir~p <-> 7r°n). For parameter set 1, we have also varied the pion 
reference energy and used either uj a or uj,. The #1,4,5,6 are insensitive to this choice, whereas 
#2,3 can vary between approximately zero and 1.5%. This points towards the necessity of a 
fourth order calculation. 





R2 [%} 


#3 [%] 


#4 [%] 


#5 [%} 


Fit 1 


0.9 


-0.5 


-0.7 


1.1 


Fit 2 


1.1 


-0.6 


-0.9 


1.1 


Fit 3 


0.9 


-0.5 


-0.8 


1.0 



Table 1: Values of the ratios Ri (i = 2, 3, 4, 5) for the various parameter sets as given by the fits 
of FMS. 

We now turn to the two isoscalar ratios. First, we consider Rq, which was first discussed by 
Weinberg [5]. For fit 1, we find Rq = 19%, which is somewhat smaller than the 25% reported 
in ref.[6]. Note, however, that the value for a + based on the KA85 phase shifts is larger (in 
magnitude) than the one used in ref.[6] (based on the LECs as determined in ref.[9]) and thus the 
isospin violating effect is indeed expected to be smaller. On the other hand, for the parameters 
of fits 2 and 3, Rq gets much larger, but not because the isospin-violating function changes 
(it is indeed stable up to fourth order as argued in ref.[6]), but rather the isospin-conserving 
function g^ 3 varies considerably. Thus, to sharpen the prediction for Rq, one has to go to next 
order in the isospin-conserving case. Interestingly, for the same parameter set (fit 1 of FMS), 
the prediction for R\ is even larger, 

R 1 = 36.7% , (20) 

which is again a huge isospin violating effect in an isoscalar quantity. The same remarks as 
made for i?i apply here. However, we stress again that this novel ratio could be accessible 
experimentally if the proposal of Bernstein [10] to extract the 7T°p scattering by precise pion 
photoproduction experiments could be carried out. Again, we remark that the value for R\ 
given in eq.(20) should be considered as a lower limit since for the other parameter sets the 
isospin-conserving isoscalar amplitude is smaller (in magnitude) which enhances the very stable 
isospin-breaking difference in the ratio. It is thus mandatory from theory and experiment to 
get a more precise value for the isoscalar amplitude. 

It is most interesting to separate the hadronic isospin violation encoded in the operator 
~ C5(m M — rrid) from the virtual photon effects. One could set C5 = 0, i.e. all strong isovector 
terms would vanish. This is not quite what one wants since then the proton is heavier than the 
neutron. In that case, the masses of the external particles would also change and a meaningful 
comparison becomes difficult. We can, however, keep C5 / for the nucleon mass insertions and 
set it to zero in the NNtttt terms ~ C5, cf. eq.(l), which we denote by ^(inr) = in table 2. 
This splitting of the strong isospin violating terms is similar to what is called "static" and 
"dynamical" isospin breaking in ref.[10]. We see that while the em isospin breaking is generally 
dominant (with the exception of Rq), there is also some sizeable strong isospin breaking. 
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Ri 


R2 


^3 


R4 


R 5 


Rq 


c 5 /0 

C 5 (7T7r) = 


36.7 
45.4 


0.9 
1.6 


-0.5 
0.8 


-0.7 
-0.7 


1.1 
1.1 


19.3 




Table 2: Values of the ratios Ri (i = 1,2,3,4,5,6) for the parameters of fit 1 from FMS 
including the full contribution from the strong isospin breaking (05 7^ 0) and with the strong 
isospin violation only contributing to the proton-neutron mass difference {c^(im) = 0). 



So far, we have discussed the ratios Ri from the real parts of the complex-valued ttN am- 
plitudes evaluated at the pertinent threshold kinematics. It is also instructive to give the 
corresponding scattering lengths. We define these as the amplitude at threshold (i.e. in 
some cases we have complex numbers) including the normalization factors, symbolically a = 
yJmxm^T^j (4iTy/s) . For scattering pions off protons, we get the numbers given in table 3. The 
imaginary part in a(ir~p — ► tt~p) is due to the intermediate ir°n state, whereas for a(-7T°p — ► -K + n) 
the energy of the initial two-body system is smaller than the one of the final two-body state. 
In both cases, these imaginary parts are fairly small. As a check on the numerics, we find that 
a(ir + n — > 7r°p) = a(ir°p — ► ir + n) and a(ir~p — ► 7r°n) = a(ir°n — ► ir~p) as demanded by time re- 
versal invariance. For comparison, the scattering lengths in the isospin limit (using the charged 
pion and the proton mass) can be obtained from table 2 of FMS by use of the relations a^^p — > 
7r ± p) = ag~ + T Oq+, a(vr°p — ► 7r°p) = »o+ an< ^ a ( 7r ~P — ^ tt ^) = a(vr°p — ► 7r + n) = — -v/2ao+- 





a(ir + p — > 7T + p) 


a(7T p — > 7T p) 


a(7T°p — > 7T°p) 


a(7r p — > 7r°n) 


a(7r°p — > 7r + n) 


Fit 1 


-108.7 


70.2 + i 3.65 


-13.4 


-125.5 


-124.7 -i 0.63 


Fit 2 


-83.8 


71.3 + i 3.65 


-0.1 


-108.6 


-107.9 - i 0.63 


Fit 3 


-94.9 


77.7 + i 3.65 


-1.8 


-121.1 


-120.2 - i 0.63 



Table 3: Values of the scattering lengths for pion scattering of protons in units of 10 3 /M 7r + for 
the various parameter sets as given by the fits of FMS. 



7. In this paper, we have considered isospin violation in low energy pion-nucleon scattering in 
the framework of heavy baryon chiral perturbation theory to third order in small momenta. We 
have taken into account all operators related to strong isospin breaking and the electromagnetic 
ones, which lead to the pion and nucleon mass differences. Stated differently, the finite parts 
of some of the virtual photon operators contributing at this order have been set to zero. This 
allows in particular to isolate the contribution of the strong dimension two isovector operator 
first considered by Weinberg. We have considered a set of six ratios Ri, which vanish in the limit 
of isospin conservation. From these, six involve isovector and isoscalar amplitudes (#2,3,4,5,7,8) 
and the two others are purely of isoscalar type (Ri t e). While in the first case, isospin violation is 
typically of the order of one percent, more sizeable effects are found in Rq [5, 6] and, as for the 
first time noted here, in R±. These results strongly motivate efforts to measure more precisely the 
isoscalar scattering length a + and try to determine the 7r°p scattering length e.g. from accurate 
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threshold pion photoproduction experiments. We also stress again that within the framework 
presented here, a unique and unambiguous separation of all different isospin violating effects is 
possible. To access the size of isospin violation encoded in the presently available pion-nucleon 
scattering data, the extension of this scheme to include Coulomb (hard) and soft photons is 
mandatory. Once this is done, it will be possible to analyze directly the cross section data 
without recourse to any model for separating em or hadronic mass effects thus, avoiding any 
mismatch by combining different approaches or models. Work along such lines is underway. 
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Figures 
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Figure 1: Tree graphs contributing to isospin violation in ttN scattering. Solid and dashed 
lines denote nucleons and pions, in order. The circle-cross (heavy dot) refers to a dimension 
two insertions ~ fa or ~ C5 (~ 1/m). 
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Figure 2: One loop graphs contributing to isospin violation in irN scattering. Solid and 
dashed lines denote nucleons and pions, in order. 
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